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Purpose  and  Conclusions 

The  purpose  of  this  report  is  to:  1)  determine  heating  value  variability  on  a  formation-by-formation 
basis;  2)  determine  the  measurement  uncertainty  due  to  gas  variability;  3)  identify  the  critical  factors 
that  cause  increased  variability;  4)  establish  gas  sampling  frequency  requirements  that  will  minimize  the 
risk  of  mis-measurement  and  mis-reporting;  and,  5)  quantify  the  potential  fiscal  impacts  of  gas 
variability.  This  report  is  not  intended  to  quantify  the  uncertainty  due  to  sampling  and  analysis 
methodology. 

The  conclusions  are  as  follows: 

•  Computer  simulation  and  the  data  collected  show  that  heating  value  variability  is  random  rather 
than  cyclic  and  shows  no  dependence  on  any  variable  analyzed  (age,  type,  formation,  pressure , 
temperature,  etc.); 

•  The  uncertainty  in  heating  value  and  relative  density  for  any  given  gas  sample  taken  does  not 
depend  on  the  frequency  of  sampling;  it  only  depends  on  the  historical  variability  of  the  heating 
value  and  relative  density; 

•  The  uncertainty  in  heating  value  and  relative  density  averaged  over  a  period  of  time  is 
dependent  on  sampling  frequency  and  an  equation  was  developed  relating  average  heating 
value  and  relative  density  uncertainty  to  sampling  frequency; 

•  Based  on  current  requirements  (i.e.  annual  sampling),  the  uncertainty  in  royalty  due  to  volume 
measurement  and  sampling  is  ±5.09%,  or  ±$127  million,  based  on  2008  revenues; 

•  There  are  numerous  options  on  how  to  address  sampling  frequency  to  minimize  the  resulting 
uncertainty  in  heating  value. 

Background 

On  December  17,  2007,  the  U.S.  Department  of  Interior,  Royalty  Policy  Committee,  Subcommittee  on 
Royalty  Management,  published  a  report  entitled  "Mineral  Revenue  Collection  from  Federal  and  Indian 
Lands  and  the  Outer  Continental  Shelf"  [RPC  report].  This  report  included  over  100  recommendations  to 
improve  the  way  Federal  and  Indian  mineral  revenues  are  collected  and  verified.  Many  of  these 
recommendations  relate  to  the  determination  and  reporting  of  gas  heating  value.  Recommendation  3-1 
states: 


“MMS  and  BLM  should  develop  a  procedure  to  determine  the  potential  BTU  variability  of  produced 
natural  gas  on  a  by-reservoir  or  by-lease  basis ,  and  estimate  the  implications  for  royalty 
payments." 

In  addition,  one  of  the  milestones  for  the  completion  of  Recommendation  3-1  is: 

"Fluids  Team  to  use  the  [variability]  report  to  develop  gas  sampling  frequency  requirements  and  to 
analyze  the  effect  of  BTU  and  gravity  variation  on  overall  measurement  uncertainty." 
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Recommendation  3-2  goes  on  to  state: 


"...MMS  and  BLM  should  determine  what  the  proper  BTU  sampling  and  analysis  frequency  should 
be  and  develop  language  for  the  revision  of  Onshore  Order  5  and  MMS  Regs.  250.1203..." 

This  report  is  an  attempt  to  address  these  recommendations. 

Before  discussing  heating  value  variability  in  any  detail,  it  is  important  to  emphasize  the  significance  of 
heating  value  on  royalties.  For  virtually  all  Indian  and  Federal  gas  produced,  royalty  due  on  production  is 
determined  with  the  following  equation: 

$R  =  RR  x  Vx  HVx  P  (Eq.  1) 


Where: 


$R  is  the  royalty  due  on  gas  produced  from  a  given  case  for  a  given  month 
RR  is  the  lease  royalty  rate  (usually  12.5%) 

V  is  the  volume  of  gas  produced  (Mcf) 

HV  is  the  heating  value  of  the  gas  that  is  produced  (Btu/scf) 

P  is  the  price,  or  monetary  value  of  the  gas  determined  by  MMS  ($/MMBtu) 

While  the  volume  (V)  produced  from  a  lease  or  case  is  determined  by  a  permanently  installed  meter 
that  must  meet  a  variety  of  standards,  the  heating  value  (HV)  is  typically  (on-shore)  determined  by 
taking  a  spot  sample  of  the  gas  and  having  it  analyzed.  The  analysis  determines  the  composition  of  the 
gas  and  then  uses  a  standard  equation1  to  calculate  the  heating  value.  In  addition  to  heating  value,  the 
gas  analysis  is  also  used  to  calculate  the  relative  density  and  supercompressibility  of  the  gas;  both  of 
which  are  used  by  flow  meters  in  the  determination  of  volume. 

While  BLM  has  standards  in  place  for  the  volume  measuring  equipment,  BLM's  only  requirement  for  gas 
sampling  is  that  it  be  done  once  per  year2.  BLM  has  no  requirements  for  sampling  methods,  sampling 
locations,  analysis  methods,  or  reporting  basis.  Referring  again  to  Equation  1,  it  can  be  seen  that 
uncertainty  in  the  heating  value  will  have  the  same  affect  on  royalty  as  uncertainty  in  volume. 

Uncertainty  in  heating  value  comes  from  two  sources:  1)  the  equipment  and  methodology  used  to 
obtain,  analyze,  and  report  heating  value,  and  2)  statistical  influences  when  taking  spot  samples.  This 
report  focuses  on  the  statistical  influences  of  spot  sampling.  Inclusion  of  sampling  and  analysis 
standards  in  the  rewrite  of  Onshore  Order  5  will  address  the  equipment  and  methodology  aspects. 

Data  Collection 

In  order  to  determine  heating  value  variability,  it  was  first  necessary  to  obtain  a  number  of  gas  analyses 
that  would  provide  the  statistical  basis  for  the  study.  Federal  Points  of  Measurement  (POM)  were 


1  API  Chapter  14.5,  “Calculation  of  Gross  Heating  Value,  Relative  Density,  Compressibility  and  Theoretical 
Hydrocarbon  Liquid  Content  for  Natural  Gas  Mixtures  for  Custody  Transfer",  API,  Third  Edition,  January,  2009 

2  Onshore  Oil  and  Gas  Order  Number  5,  1989,  III.C.23 
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chosen  that  provided  a  reasonable  range  of  reservoir  types,  flowing  conditions,  surface  and  downhole 
plumbing  configurations,  and  age.  In  May,  2009,  Written  Orders  to  the  operators  of  the  selected  POMs 
were  sent  out  asking  for  the  gas  analyses  for  the  past  6  years  or  the  life  of  the  meter,  whichever  was 
longest.  Table  1  summarizes  the  records  requested  and  the  records  received. 


Field  Office 

#  of  POM 
requested 

#  of  POM 
received 

#  of  POM 
used 

#  of  samples 

Buffalo,  WY 

62 

42 

42 

221 

Farmington,  NM 

65 

64 

59 

549 

Glenwood  Springs,  CO 

24 

9 

9 

24 

Pinedale,  WY 

81 

52 

46 

513 

Rawlins,  WY 

46 

38 

37 

292 

Vernal,  UT 

44 

32 

24 

296 

TOTAL 

322 

237 

217 

1895 

Table  1  -  Data  summary 

As  shown  in  Table  1,  certain  records  were  received  that  could  not  be  tied  to  a  known  lease  or  case 
number,  or  were  not  included  in  the  original  request.  Because  the  conditions  of  these  meters  were  not 
known,  these  records  could  not  be  used  for  most  of  the  analysis. 

Once  the  data  were  received,  they  were  key  entered  into  an  Excel  spreadsheet  and  then  imported  into 
an  Access  database  for  manipulation  and  analysis.  Some  data  filtering  and  editing  was  also  performed. 

Heating  Value  Variability  By  Formation 

This  analysis  is  intended  to  address  part  of  Recommendation  3-1  in  the  RPC  report,  which  asks  BLM  to 
develop  a  procedure  to  determine  BTU  variability  of  produced  natural  gas  on  a  by-reservoir  or  by-lease 
basis.  It  was  determined  that  there  was  insufficient  data  to  perform  an  analysis  on  a  lease-by-lease 
basis.  BLM  currently  manages  about  23,000  producible  oil  and  gas  leases.  The  data  for  237  POMs 
received  (see  Table  1)  covers  193  different  leases.  Therefore,  the  data  for  this  analysis  represents  less 
than  1%  of  the  total  number  of  leases  that  BLM  manages,  making  it  virtually  impossible  to  perform  an 
analysis  on  a  lease-by-lease  basis. 

It  was  also  not  possible  to  perform  the  analysis  on  a  reservoir-by-reservoir  basis  because  of  the  difficulty 
in  defining  a  "reservoir"  and  in  obtaining  data  by  reservoir.  Instead,  it  was  decided  to  perform  the 
analysis  on  a  formation-by-formation  basis  as  this  is  a  better  defined  parameter  and  the  producing 
formation  data  are  readily  available  to  BLM. 

The  gas  analysis  data  base  was  first  sorted  on  formation,  meter  name,  and  sample  date.  For  each 
formation,  the  average  heating  value  of  all  the  samples  taken  from  that  formation  was  calculated  along 
with  the  standard  deviation.  The  scatter  in  the  heating  values  was  determined  at  the  95%  confidence 
level  using  a  student  t  -distribution.  For  example,  the  Almond  formation  has  an  average  heating  value  of 
1109.5  Btu/scf  with  a  95%  confidence  variation  of  ±4.91%.  This  indicates  that  95%  of  the  heating  values 
obtained  were  between  1055.1  Btu/scf  and  1164.0  Btu/scf.  Use  of  the  student  t-distribution  allows  for 
meaningful  analysis  of  small  data  sets.  Table  2  summarizes  the  this  analysis  (note  that  variation  is  not 
calculated  when  there  are  less  than  5  samples). 
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Average 


Heating 

Variation 

Formation 

#meters  #samples 

Value 

95% 

(Btu/scf) 

(±%) 

ALMOND 

6 

31 

1109.55 

4.91 

ALMY 

4 

29 

1186.04 

11.03 

ANDERSON  COAL 

5 

28 

902.65 

19.12 

BAXTER 

3 

15 

1115.99 

9.40 

BEAR  RIVER 

4 

4 

1113.62 

BEAR  RIVER/FRONTIER  2 

1 

52 

1126.95 

2.64 

BIG  GEORGE  COAL 

5 

21 

945.39 

1.99 

CANYON  COAL 

3 

18 

942.64 

15.48 

CANYON  COAL  (LOWER) 

1 

1 

981.60 

CARMEL 

1 

11 

1037.74 

1.13 

CASTLEGATE 

2 

24 

1141.63 

9.07 

CEDAR  MOUNTAIN 

1 

3 

1029.75 

CHACRA 

5 

31 

1208.99 

11.86 

COOK COAL 

2 

6 

908.95 

48.30 

COW  CREEK 

1 

5 

1118.44 

5.57 

DAKOTA 

9 

47 

1240.95 

17.48 

DAKOTA /MORRISON 

1 

7 

1039.01 

0.94 

DEEP  CREEK  1  1 

988.49 

DIETZ  3  COAL 

1 

4 

957.71 

ERICSON 

1 

4 

1115.10 

FARMINGTON 

2 

10 

1056.70 

8.67 

FELIX  COAL 

1 

4 

920.44 

FORT  UNION 

3 

22 

923.30 

5.32 

FORT  UN  ION  COAL 

2 

25 

958.77 

3.96 

FRONTIER 

11 

78 

1093.39 

7.53 

FRONTIER  2 

2 

9 

1123.27 

4.72 

FRUIT  LAND 

9 

104 

1044.44 

15.08 

FRUIT  LAND  COAL 

1 

22 

956.29 

1.14 

GALLUP 

4 

27 

1080.17 

17.75 

GREEN  RIVER 

2 

11 

1041.22 

8.93 

HILLIARD/BAXTER 

1 

47 

1160.95 

3.65 

KIRTLAND 

1 

2 

967.20 

KNOBLOCH  COAL 

1 

18 

938.69 

0.94 

LANCE 

2 

8 

1190.49 

323 

LANCE  POOL 

5 

29 

1083.96 

1.73 

LANCE  POOLTAN SI  LACE  P®L 

1 

13 

1089.82 

0.43 

LANCE/WARDELL 

1 

11 

1145.56 

0.77 

LEWIS 

11 

87 

1143.54 

9.33 

LOWER  CANYON  COAL 

1 

10 

944.11 

2.14 

LOWER  COOK  COAL 

1 

4 

971.58 

MAN COS  B 

1 

4 

1314.96 

MESAVERDE 

26 

213 

1146.28 

9.06 

MESAVERDE 1 

5 

44 

1184.63 

19.28 
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Average 

Heating  Variation 

Formation  #meters  #samples  Value  95% 

(Btu/scf)  (±%) 


MESAVERDE  1/LANCE 

1 

6 

1085.00 

0.57 

MIDWEST  COAL 

2 

2 

893.77 

MONARCH  COAL 

1 

3 

979.67 

MORGAN 

1 

4 

993.49 

MORRISON 

1 

12 

1110.63 

0.73 

MUDDY 

3 

18 

1155.64 

24.60 

NACIMIENTO 

1 

38 

862.14 

5.53 

NIOBRARA 

1 

2 

1420.13 

OJO  ALAMO 

3 

33 

1041.75 

2.06 

PARKMAN 

1 

3 

1593.47 

PICTURED  CLIFFS 

7 

71 

1154.36 

15.47 

PINON  FRUITLAND  SAND 

1 

5 

1171.40 

1.35 

POINT  LOOKOUT 

5 

37 

1047.90 

16.05 

ROBERTSON  COAL 

1 

27 

1050.24 

13.83 

SI-FRONTIER  1/S2-BEAR  RIVER 

1 

7 

1150.09 

1.71 

SI-FRONTIER  1/S2-FRONTIER  2 

2 

13 

1105.29 

1.57 

S1-FRONTIER-MUDDY/S3-BAXTER 

1 

5 

1140.63 

2.01 

SAN  JOSE 

3 

26 

1108.75 

12.85 

SHANNON 

2 

8 

1163.82 

13.70 

SKULL CREEK 

1 

3 

1300.97 

SMITH  COAL 

1 

2 

987.00 

TANSILL 

1 

5 

1179.81 

1.35 

TOCITO 

1 

19 

1347.50 

18.73 

TRANSITION 

2 

26 

1242.82 

4.69 

TRANSIT  ION/MESAVERDE  1 

1 

5 

1198.03 

9.68 

TURNER 

3 

15 

1387.71 

8.69 

UINTA 

1 

11 

1015.09 

0.12 

UPPER  COOK  COAL 

1 

6 

953.22 

9.72 

WALL  COAL 

1 

2 

935.59 

WASATCH 

9 

75 

1142.01 

17.76 

WASATCH G SAND 

1 

2 

1069.87 

WILLIAMS  FORK 

2 

11 

1118.92 

6.36 

WOLF  CAMP 

1 

4 

1113.58 

WYODAKCOAL 

1 

12 

901.45 

6.34 

Table  2  -  Heating  Value  Variation  by  Formation 


Sampling  Methods 

Referring  to  Equation  1,  royalty  payments  are  directly  affected  by  the  measured  heating  value  of  a  gas. 
A  10%  error  in  the  measured  heating  value  will  result  in  a  10%  error  in  royalty  payment  as  long  as  the 
errant  gas  analysis  report  is  in  effect. 
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It  is  not  unreasonable  to  expect  that  heating  value  will  vary  over  time  depending  on  geographic  location 
within  a  formation,  changing  reservoir  conditions,  changing  ambient  conditions,  changing  flowing 
conditions,  or  changing  gas  streams  for  meters  that  measure  more  than  one  well.  Stratification  or 
segregation  of  heavier  and  lighter  gas  can  occur  in  a  single  formation  causing  variation  in  the  gas  stream 
over  time  or  as  a  function  of  producing  depth  or  location  of  the  well  within  the  formation. 

To  account  for  changes  in  heating  value  over  time,  only  three  different  measurement  methodologies 
have  been  adopted  by  industry:  on-line  gas  chromatographs,  composite  sampling,  and  spot  sampling 
(cylinder  and  portable  gas  chromatograph). 

One  method  of  determining  heating  value  is  through  continuous  measurement,  which  can  be 
approximated  with  an  on-line  gas  chromatograph.  This  device  automatically  samples  and  analyzes  the 
gas  stream  every  few  minutes  and  then  transmits  the  calculated  heating  value,  relative  density,  and 
composition  to  the  flow  computer.  The  flow  computer  multiplies  each  measured  heating  value  by  the 
volume  of  gas  measured  to  integrate  the  energy  content  of  the  gas  over  a  reporting  period.  On  line  gas 
chromatographs  are,  however,  expensive  and  are  typically  only  found  on  high-volume  sales  meters. 

A  second  method  of  determining  heating  value  is  with  a  flow-proportional  composite  sample.  This  is  a 
device  that  allows  a  small  sample  of  gas  to  be  added  to  a  sample  cylinder  at  even  increments  of  volume 
measured  by  the  flow  meter.  Once  the  sample  cylinder  is  filled,  it  is  taken  to  a  lab  where  the  gas  is  put 
through  a  gas  chromatograph  to  determine  the  composition.  Over  a  period  of  time,  the  sample  cylinder 
will  contain  a  true  volume-weighted  average  composition  from  which  heating  value  and  relative  density 
are  calculated.  As  with  on-line  gas  chromatographs,  these  devices  also  tend  to  be  found  only  on  higher 
volume  sales  meters. 

The  third  form  of  sampling  is  with  a  spot  sample  taken  at  some  interval  of  time  (minimum  of  once  per 
year  for  BLM).  A  spot  sample  is  taken  by  either  filling  a  sample  cylinder  and  taking  it  to  a  lab  for  analysis 
or  bringing  a  portable  gas  chromatograph  to  the  field  and  temporarily  connecting  it  directly  to  the 
sample  port.  Through  inquiries  with  BLM  Field  Offices,  it  is  estimated  that  at  least  95%  of  heating  values 
for  onshore  Federal  and  Indian  gas  are  determined  through  spot  sampling. 

Uncertainty  and  Bias  Due  to  Spot  Sampling 

The  uncertainty  of  heating  value  and  relative  density  is  dependent  on  the  sample  and  analysis 
equipment  and  methodology  for  all  sampling  methods  described  above.  Sample  bias  is  minimized  by 
following  published  industry  standards  that  are  based  on  research  and  field  experience3.  However, 
because  spot  sampling  is  essentially  a  random  snapshot  in  time,  it  also  includes  a  statistical  component 
of  uncertainty.  This  "statistical"  source  of  uncertainty  is  the  focus  of  this  study.  Total  uncertainty  in 
heating  value  (UHV)  taken  by  spot  sampling  is  the  combination  of  the  uncertainty  due  to  the  equipment 
and  procedures  used  (Uep)  and  the  statistical  uncertainty  (Us).  In  equation  form,  total  uncertainty  is 
determined  as  follows: 


(Eq.  2) 


3  Although  BLM  does  not  currently  incorporate  or  enforce  these  standards,  this  is  beyond  the  scope  of  this  report 
and  is  being  addressed  in  other  RPC  recommendations. 


Online  gas  chromatographs  and  composite  samplers  effectively  integrate  the  composition  overtime  and 
do  not  include  the  statistical  component. 

To  demonstrate  the  uncertainty  caused  by  spot  sampling,  assume  that  the  true  heating  value  at  a 
particular  meter  varies  with  time  as  shown  in  Figure  1.  Also  assume  that  a  single  spot  sample  is  taken  at 
some  point  during  the  year  shown  and  the  heating  value  derived  from  the  spot  sample  is  used  for 
royalty  determination  for  the  entire  year.  From  Figure  1,  the  heating  value  used  for  royalty 
determination  could  be  anywhere  between  950  Btu/scf  and  1050  Btu/scf,  depending  on  the  time  and 
date  the  sample  was  taken.  This  gives  an  uncertainty  of  +50  Btu/scf,  or  ±5%.  This  statistical  uncertainty 
would  be  in  addition  to  the  uncertainty  due  to  the  sampling  equipment  and  procedures  used  to  obtain 
the  sample  and  analysis. 


Effect  of  Variability  on  Spot  Sampling 


Variability  patterns 


Figure  1 


One  of  the  first  hurdles  in  this  analysis  is  determining  the  type  of  variation  that  is  present.  Is  the 
variation  highly  random  and  variable  such  as  that  shown  in  Figure  1,  or  is  it  more  cyclic  as  shown  in 
Figure  2  (reflecting  a  possible  seasonal  fluctuation,  for  example).  From  Figure  2,  it  can  be  seen  that  if 
samples  were  taken  frequently  -  once  a  week  for  example  -  the  change  in  heating  value  between 
samples  would  be  minimal.  If,  on  the  other  hand,  samples  were  only  taken  once  or  twice  per  year,  the 
odds  of  a  much  larger  fluctuation  between  heating  values  would  be  greater.  To  demonstrate  this  effect, 
a  computer  simulation  was  developed  that  generated  random  cyclic  patterns  similar  to  that  shown  in 
Figure  2.  Then,  virtual  samples  were  taken  at  different  frequencies  and  the  average  difference  in  heating 
values  from  sample  to  sample  were  plotted  for  each  frequency.  The  results  of  this  simulation  are  shown 
in  Figure  3,  where  each  dot  represents  the  difference  in  heating  value  between  the  sample  taken  and 
the  previous  sample  taken.  There  is  a  distinct  relationship  between  sample  frequency  and  maximum 
change  in  heating  value  between  samples  (shown  by  the  red  line).  However,  once  the  time  between 
samples  exceeds  60  days,  or  about  a  2/3  of  the  cycle  frequency  (from  Figure  2),  the  relationship  no 
longer  applies  and  becomes  random. 
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Example  of  Cyclic  Variation 


Detection  of  Cyclic  Pattern 


Time,  days 

Figure  2 


Time  between  samples,  days 
Figure  3 


The  same  type  of  analysis  was  done  on  the  data 
collected  in  an  attempt  to  indentify  some  type  of  relationship  between  sample  frequency  and  change  in 
heating  value  between  samples.  If  a  relationship  similar  to  that  shown  in  Figure  3  is  apparent,  it  could 
indicate  that  heating  value  variation  is  cyclic  in  nature.  This  analysis  was  done  on  a  formation-by- 
formation  basis.  Graphs  of  the  four  formations  with  the  greatest  amount  of  data  are  shown  in  Figures  4a 
through  4d.  None  of  these  figures  demonstrate  any  relationship  between  sampling  frequency  and 
change  in  heating  value,  nor  did  any  of  the  other  formations  analyzed.  This  leads  to  the  conclusion  that 
there  is  no  evidence  indicating  that  heating  value  variation  as  a  function  of  time  is  cyclic  in  nature. 


It  is  possible  that  the  uncertainty  in  the  sampling  and  analysis  methods  mask  any  actual  variation  in 
heating  value  or  it  is  possible  that  the  cyclic  pattern,  if  it  exists,  is  shorter  than  the  sampling  frequencies 
analyzed.  Regardless  of  the  cause,  however,  it  is  assumed  for  the  remainder  of  the  study  that  heating 
value  variation  is  completely  random,  as  shown  in  Figure  1. 


Uncertainty  of  single  heating  values 

Going  back  to  the  pattern  in  Figure  1,  now  assume  that  instead  of  a  single  annual  sample,  a  monthly 
sample  is  taken  sometime  within  each  month  and  used  to  report  the  heating  value  on  OGOR  B4  for  that 
month.  The  range  of  values  that  could  be  obtained  for  each  sample  is  shown  by  Table  3.  From  Table  3,  it 
can  be  seen  that  increasing  sample  frequency  from  once  per  year  to  once  per  month  has  very  little 
effect  on  the  uncertainty  of  a  sample. 


To  test  this  hypothesis  -  that  sample  frequency  has  little  impact  on  individual  values  of  heating  value 
uncertainty -the  data  collected  were  again  analyzed.  First,  the  data  were  sorted  to  exclude  meters 
where  there  were  less  than  5  samples.  Then,  the  variability  of  the  heating  values  for  each  meter  was 
calculated  using  a  student  t-distribution  as  shown  by  Equation  3  (See  Appendix  1  for  an  example  of  this 
method). 


where: 


(Eq.3) 


4"OGOR  B",  is  the  Oil  and  gas  Operator  Report,  Part  B,  which  is  an  MMS  form  used  to  collect  volume  and  heating 
value  on  which  royalty  is  calculated  for  a  lease,  participating  area,  or  communitization  area. 


UHVis  the  uncertainty  (95%  confidence  level)  of  the  heating  value,  ±% 

o  is  the  standard  deviation  of  the  previous  samples,  ±% 

tdis,  is  the  "student  t-distribution"  function  shown  in  Appendix  2 


Month 

Heating  Value  Range  (Btu/scf) 

Uncertainty 

High 

Low 

(±%) 

January 

1047 

953 

4.70 

February 

1048 

953 

4.75 

March 

1050 

951 

4.95 

April 

1049 

952 

4.85 

May 

1050 

952 

4.90 

June 

1046 

958 

4.40 

July 

1049 

951 

4.90 

August 

1040 

952 

4.40 

September 

1049 

958 

4.55 

October 

1045 

957 

4.40 

November 

1047 

956 

4.55 

December 

1050 

951 

4.95 

Table  3 


Note  that  Equation  3  is  used  to  determine  the  total  uncertainty  in  heating  value  (Uhv)  because  the 
variability  in  the  historical  data  include  the  effects  of  both  statistical  uncertainty  (Us)  and  uncertainty 
due  to  the  equipment  and  procedures  used  (Uep).  Finally,  the  calculated  variability  of  the  heating  value 
for  each  meter,  using  Equation  3,  was  plotted  against  the  average  time  between  samples  for  that  meter 
as  shown  in  Figure  5.  Again,  there  is  no  apparent  dependence  of  variability  on  sampling  frequency, 
which  not  only  supports  the  hypothesis  described  above,  but  also  reaffirms  the  apparent  random  nature 
of  heating  value  variations. 

Uncertainty  of  average  heating  value 

Another  way  to  look  at  statistical  uncertainty  due  to  spot  sampling  is  the  uncertainty  in  average  heating 
value  taken  over  some  period  of  time,  rather  than  the  uncertainty  of  any  given  sample.  For  the  pattern 
shown  in  Figure  1,  a  time  was  randomly  picked  within  each  month  when  a  virtual  sample  was  taken.  This 
is  shown  by  the  yellow  dots  in  Figure  6.  Table  4  lists  the  heating  values  obtained  from  the  12  random 
monthly  samples.  From  Table  4,  the  average  heating  value  for  the  year  is  995.8  Btu/scf,  whereas  the 
true  average  of  the  data  set  is  999.9,  for  an  error  of  only  0.41%.  By  taking  12  monthly  samples,  the 
uncertainty  in  heating  value  due  to  spot  sampling  over  a  year's  time  has  been  reduced  from  ±5%  to  only 
±0.41%. 

Of  course,  had  another  pattern  been  used  in  Figures  1  and  6,  or  other  sampling  points  selected,  this 
averaging  error  could  have  been  lower  than  the  true  average,  higher  than  the  true  average,  or  even  the 
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Sample  Date 

s-6 - 

Heating  Value 
(Btu/scf) 

January  16 

971.6 

February  27 

1042.3 

March  12 

990.5 

April  29 

1016.0 

May  16 

965.0 

June  5 

1001.5 

July  12 

993.0 

August  29 

1014.0 

September  9 

1007.7 

October  16 

983.5 

November  22 

998.6 

December  16 

965.6 

Average 

995.8 

Table  4 


same  as  the  true  average.  Therefore,  this  averaging  error  actually  becomes  a  range  of  potential  error 
based  on  the  specific  case.  A  "range  of  potential  error "  is  the  definition  of  uncertainty;  therefore,  there 
is  uncertainty  associated  with  the  statistical  probability  of  spot  sampling. 

A  general  conclusion  can  be  drawn  that  the  more  samples  that  are  taken  over  a  period  of  time,  the 
closer  the  average  of  the  samples  will  be  to  the  true  average.  The  longer  the  time  between  samples,  the 
fewer  samples  are  collected.  While  the  monthly  sampling  results  in  12  data  points  for  the  year,  there 
would  only  be  1  data  point  if  an  annual  sample  was  taken.  The  fewer  data  points  there  are  in  an 
average,  the  higher  the  probability  that  the  average  will  not  accurately  represent  the  data  set.  This  is  the 
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basis  of  the  remainder  of  this  study  -  establishing  a  relationship  between  sample  frequency  and  the 
average  heating  value  uncertainty  over  a  period  of  time. 

The  selection  of  a  time  frame  for  which  an  average  is  calculated  is  somewhat  arbitrary.  A  period  of  1 
year  was  used  in  this  study.  The  selection  of  a  shorter  averaging  time  will  result  in  a  higher  uncertainty 
for  a  given  sampling  frequency  and  a  longer  averaging  time  will  result  in  lower  uncertainty  than  those 
given. 

Computer  Modeling 

To  quantify  the  uncertainty  of  average  heating  value  due  to  variability,  a  computer  "Monte  Carlo"  model 
was  developed.  The  model  first  creates  a  random  pattern  of  heating  values  over  a  1-year  span  of  time  in 
one-hour  increments,  similar  to  Figure  1.  Then,  virtual  samples  of  heating  value  are  taken  from  the 
pattern  at  sampling  intervals  ranging  from  once  per  day  to  once  every  360  days,  in  one  day  steps.  An 
average  heating  value  for  the  pattern  is  determined  and  recorded  for  each  sampling  frequency.  Once 
that  is  complete,  a  new  pattern  is  generated  and  samples  are  once  again  taken.  This  process  is  repeated 
100  times  to  determine  the  relationship  between  sampling  frequency  and  uncertainty. 

For  each  run  of  the  model,  the  maximum  amount  of  variation  along  with  the  central  value  can  be 
assigned.  A  central  value  of  1000  Btu/scf  was  used  for  all  runs,  and  four  levels  of  maximum  variation 
were  used  in  the  analysis:  ±500  Btu/scf,  ±100  Btu/scf,  ±50  Btu/scf,  and  ±10  Btu/scf.  A  variation  of  ±0 
Btu/scf  was  also  run  to  ensure  the  average  and  standard  deviation  calculations  were  performed 
correctly. 

Model  results  include  the  average  heating  value 
for  each  sampling  frequency  for  each  run  of  the 
model.  Figure  7  shows  the  results  from  10  runs 
(100  runs  was  too  cumbersome  to  plot)  for  a 
maximum  assigned  pattern  variation  of  ±50 
Btu/scf.  From  Figure  7,  it  can  be  seen  that  for 
very  short  sampling  frequencies,  the  average 
heating  value  taken  over  a  year  has  much  less 
scatter  than  for  samples  of  the  same  data  taken 
less  frequently.  For  example,  the  range  of 
possible  average  heating  values  for  a  30-day 
sampling  frequency  is  about  ±8  Btu/scf, 
whereas  the  range  is  about  ±30  Btu/scf  when 
samples  are  taken  annually.  Clearly,  there  is  a 
relationship  between  sample  frequency  and 
range  of  error,  or  uncertainty,  of  average 
heating  value  taken  over  a  span  of  time.  This  is 
a  statistical  uncertainty  that  is  unique  to  spot  sampling  as  it  has  nothing  to  do  with  the  methods  or 
equipment  used  to  obtain  heating  values  in  the  field. 

To  quantify  this  relationship,  the  results  of  the  model  were  sorted  by  sampling  frequency  and  by  run  as 
shown  in  Table  5.  For  each  row  of  data,  the  amount  of  scatter  for  that  sampling  frequency  was 
calculated  as  two  times  the  standard  deviation  of  the  values  (2o).  For  a  large  data  set,  a  2o  value 
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Sample 

Frequency 

(days) 

Average  Pleating 

Value  (Btu/scf)  for  Model  Run  #... 

2a 

Variation 

(±%) 

1 

2 

3 

4 

5 

6 

7... 

100 

1 

1000.1 

1000.3 

1000.0 

999.7 

1000.2 

999.9 

999.6 

1000.1 

0.1615 

2 

999.5 

999.7 

1000.2 

100.0 

999.8 

1000.0 

999.4 

999.9 

0.2284 

3... 

999.9 

1000.0 

999.7 

1000.9 

1000.2 

999.4 

999.4 

1000.1 

0.2797 

360 

973.5 

984.0 

1010.8 

1001.5 

1024.9 

998.7 

1003.3 

989.0 

3.0642 

Table  5-  Sample  Model  Results  for  a  Random  Variation  of  ±50  Btu/scf 


yields  a  95%  confidence  level;  in  other  words,  95%  of  the  data  is  between  the  (+)  and  (-)  value  given. 
Similar  calculations  were  done  with  the  random  pattern  variation  set  at  ±100  Btu/scf  and  ±10  Btu/scf. 
Predictably,  the  higher  the  variation  set  in  the  model,  the  greater  the  2 a  variation  for  a  given  sample 
frequency  and  visa  versa.  Figure  8  shows  a  graph  of  2a  variation  (uncertainty)  as  a  function  of  sample 
frequency  and  degree  of  random  pattern  variation.  The  blue  dots  represent  the  2o  variation  in  the 
model  output  as  shown  in  Table  5  and  the  red  lines  are  curve  fits  of  each  set  of  data.  The  2o  variation  of 
the  actual  random  heating  value  patterns  was  also  calculated  (V95*). 


Sampling  Uncertainty,  2 a 


Sampling  Period  (Ps),  days 

Figure  8 


A  general  equation  (shown  by  the  red  lines  in  Figure  8)  for  each  data  set  was  determined  as  follows: 


(Eq.4) 


Where: 


Us  is  the  statistical  spot  sampling  uncertainty 

C  is  a  constant  that  is  a  function  of  the  data  variability,  V95%  (see  Table  6) 
Ps  is  the  time  between  samples,  days 
Pa  is  the  averaging  period,  days 
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Variation  in 
Random  Pattern 
(±Btu/scf) 

2o  Variation  in 

Random  Pattern 

(±%) 

"C"  for 
Equation  2 

100 

11.54 

1.097 

50 

5.772 

5.489 

10 

1.154 

10.975 

Table  6 


It  was  then  determined  that  the  constants,  "C",  from  Table  6  had  a  linear  relationship  with  the  2o 
variation  in  the  random  patterns  generated  by  the  model.  Therefore,  a  single  equation  relating 
statistical  sampling  uncertainty  (Us)  to  sampling  period  (Ps),  averaging  period  (Pa),  and  2o  variation  in  the 
sampling  data  (V95*)  was  derived  as  shown  by  Equation  5: 


(Eq.  5) 


Application  of  the  Computer  Model 

Use  of  the  computer  model  demonstrates  that  uncertainty  in  average  heating  value  due  to  sampling 
frequency  is  predictable  and  quantifiable.  To  apply  the  uncertainty  equation  (Eq.  5),  the  primary 
unknown  is  the  variability  of  the  heating  value  data  (V95*),  which  can  be  determined  using  historical 
heating  values  and  standard  statistical  methods  described  by  Equation  3  and  as  shown  in  Appendix  1. 
Once  this  has  been  determined,  Equation  5  can  be  used  to  determine  the  uncertainty  in  average  heating 
value  at  various  sampling  frequencies.  For  a  variation  (V95%)  of  ±2.47%  from  the  example  in  Appendix  1, 
the  sampling  uncertainty  at  various  sampling  frequencies  can  then  be  calculated  as  shown  in  Table  7 
(assume  a  1-year  averaging  period). 

Because  historical  data  is  used  to  determine  the  variability,  and  the  historical  data  includes  the 
uncertainty  effects  of  both  statistical  variability  (Us)  and  the  equipment  and  procedures  used  (Uep),  the 
sampling  uncertainty  obtained  using  Equation  5  is  the  total  uncertainty  in  heating  value  (U„v). 


Sample  Frequency 
(days) 

V9S* 

1-year  average  Fleating 
Value  Uncertainty 
(±%) 

30 

0.2867 

2.47 

0.668 

60 

0.4054 

2.47 

0.945 

90 

0.4966 

2.47 

1.157 

180 

0.7022 

2.47 

1.636 

365 

1.0000 

2.47 

2.330 

Table  7 


This  methodology  was  applied  for  each  meter  measuring  gas  from  the  formations  presented  in  Table  1. 
Appendix  3  shows  the  heating  value  uncertainties  of  each  meter  as  a  function  of  sampling  frequency. 

The  exact  same  methodology  can  be  used  to  determine  the  uncertainty  in  relative  density  as  follows: 

(Eq.  6) 
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where  UGris  the  uncertainty  in  relative  density.  Appendix  4  shows  the  relative  density  uncertainties. 
Meters  where  there  are  less  than  five  samples  were  excluded  from  the  analyses. 

Royalty  Implications 

Referring  to  Equation  1,  royalty  is  based  on  both  volume  and  heating  value.  Errors  in  either  value  will 
result  in  proportional  errors  in  royalty.  Likewise,  uncertainty  in  volume  and  heating  value  will  result  in 
uncertainty  in  royalty.  For  the  purpose  of  this  study,  it  will  be  assumed  that  biases  in  both  volume  and 
heating  value  have  been  minimized  by  following  industry  standards.  This  assumption,  however,  is  not 
currently  applicable  because  BLM  does  not  have  any  standards  for  sampling  and  analysis. 

While  the  focus  of  this  study  has  been  on  heating  value  uncertainty,  this  analysis  and  uncertainty 
determination  method  applies  equally  to  relative  density  determination,  as  shown  in  Appendix  4. 
Relative  density  is  important  because  it  is  used  in  the  determination  of  volume.  Therefore,  uncertainty 
in  sampling  results  in  uncertainty  to  both  volume  and  heating  value. 

In  addition  to  relative  density,  sampling  uncertainty  also  affects  the  calculation  of  supercompressibility, 
which  also  affects  volume  calculation.  However,  because  the  uncertainty  effects  on  supercompressibility 
are  minor  compared  with  relative  density,  this  effect  is  not  included  in  this  study. 

Based  on  the  data  and  results  in  Appendices  3  and  4,  an  average  uncertainty  in  both  heating  value  and 
relative  density,  based  on  the  current  requirement  for  annual  sampling,  was  determined  as  follows: 

Uncertainty  in  heating  value  (UHV):  ±3.58% 

Uncertainty  in  relative  density  (UGr):  ±4.03% 

In  addition,  the  uncertainty  in  volume,  excluding  the  effects  of  relative  density,  is  currently  limited  to 
±3%  for  meters  measuring  more  than  100  Mcf/day,  through  the  various  statewide  Notices  to  Lessees 
(NTLs).  Overall  uncertainty  of  both  volume  and  heating  value,  which  would  also  reflect  the  uncertainty 
in  royalty  is  the  root-sum-square  of  these  three  sources  of  uncertainty  and  is  determined  as  follows: 


Where: 


Uf  is  the  uncertainty  in  royalty,  ±% 

Xv  is  the  sensitivity  coefficient  for  volume  uncertainty  (1.0) 
Uv  is  the  uncertainty  in  volume,  ±% 

Xd  is  the  sensitivity  coefficient  for  relative  density  (0.5)5 
Ud  is  the  uncertainty  in  relative  density,  ±% 

XHV  is  the  sensitivity  coefficient  for  heating  value  (1.0) 

UHV  is  the  uncertainty  in  heating  value,  ±% 


=  ±5.09% 


5  In  the  flow  equation,  the  relative  density  is  taken  to  the  0.5  power;  sensitivity  coefficient  is  generally  taken  as  the 
exponent  of  the  variable. 
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Given  the  2008  total  onshore  Federal  and  Indian  royalty  income  of  about  $2.5  billion,  the  total 
uncertainty  due  to  volume  and  heating  value  is  ±$127  million,  based  on  a  1-year  average. 

If,  for  example,  the  sampling  frequency  was  increased  to  monthly  for  all  meters  used  in  this  study,  the 
uncertainty  in  royalty  would  be  reduced  from  ±5.09%  to  ±3.22%,  thereby  reducing  the  uncertainty  in 
royalty  to  ±$81  million.  Going  to  composite  sampling  or  on-line  gas  chromatographs  could  reduce 
sampling  and  analysis  uncertainty  further. 

Establishing  Minimum  Sampling  Frequency 

In  the  development  of  Onshore  Order  5,  it  will  be  necessary  to  establish  sampling  frequencies  and 
methods  that  minimize  the  risk  of  mis-measurement,  or  uncertainty,  of  heating  value  and  relative 
density.  Onshore  Order  5,  III.C.23  currently  requires  annual  sampling  for  all  sales  and  allocation  meters. 
From  Appendices  3  and  4,  it  can  be  seen  that  annual  sampling  leads  to  high  levels  of  uncertainty  in  many 
situations. 

There  are  a  number  of  options  to  consider  as  described  below: 

1.  Establish  prescriptive  sampling  frequencies  based  on  volume.  For  example,  meters  measuring 
more  than  100  Mcf/day  could  be  sampled  quarterly,  whereas  meters  measuring  100  Mcf/day  or 
less  could  be  sampled  annually.  While  this  option  would  be  easy  to  implement,  it  does  not 
necessarily  achieve  the  desired  results  of  minimizing  uncertainty.  For  example,  from  Appendix  3, 
a  quarterly  sampling  could  produce  very  low  uncertainty  for  well  "Palm  Tree  Fed  34-4",  in  the 
Big  George  Coal  (±0.26%),  but  would  result  in  a  relative  high  uncertainty  for  well  "20-mile  23A- 
1313",  in  the  Anderson  Coal  (±14.83%). 

2.  Establish  sampling  frequencies  based  on  volume  and  formation.  For  example,  meters  measuring 
gas  from  the  Anderson  Coal,  over  a  certain  volume  threshold,  would  have  to  sample  monthly, 
whereas  meters  measuring  gas  from  the  Big  George  Coal  would  only  be  required  to  sample 
annually.  This  method  would  be  effective  and  relative  easy  to  enforce.  Unfortunately,  it  would 
also  take  a  great  deal  more  research  as  the  values  in  Table  5  may  not  represent  the  true 
variability  of  each  formation. 

3.  Establish  an  uncertainty  limit  for  heating  value,  similar  to  the  uncertainty  limit  already 
established  for  volume  measurement.  While  this  approach  would  ensure  that  BLM  regulates 
sampling  uncertainty,  it  could  be  difficult  to  enforce.  For  each  point  of  measurement,  the 
historical  heating  values  would  be  needed  to  implement  the  methodologies  presented  in  this 
study.  Volume  measurement  uncertainty,  by  contrast,  can  be  determined  with  a  single  field 
visit.  Implementing  this  option  would  require  that  all  samples  from  all  Federal  and  Indian  points 
of  measurement  be  submitted  to  BLM  and  entered  into  a  database.  With  this  in  place,  a  BLM 
inspector  could  easily  find  the  heating  value  and  relative  density  uncertainty  for  a  specific 
meter. 

4.  Establish  an  uncertainty  limit  for  energy  rather  than  having  individual  uncertainties  for  volume 
and  heating  value.  This  approach  would  be  the  most  effective  method  because  royalties  are 
ultimately  based  on  energy  and  not  just  volume  or  heating  value.  It  would  also  offer  the  most 
flexibility  for  industry  as  they  could  make  up  for  a  high  sampling  uncertainty  with  lower  volume 
uncertainties.  Flowever,  this  would  require  a  data  base  explained  in  the  number  3,  above. 
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5.  Establish  an  initial  sampling  frequency  that  could  later  be  modified  once  there  were  enough  data 
to  perform  an  analysis.  The  initial  frequency  could  be  relatively  restrictive  -  once  a  month,  for 
example,  and  put  the  burden  on  industry  to  provide  the  data  to  demonstrate  a  less  frequent 
sampling  is  warranted,  based  on  the  methodology  presented  in  this  study.  This  would  be  on  a 
case-by-case  or  area-by-area  basis.  Downsides  to  this  approach  would  be  the  creation  of  a  large 
workload  for  BLM  in  reviewing  the  data  submitted  and  approving  and  tracking  lesser  sampling 
frequencies. 


The  BLM  Onshore  Order  5  team  will  be  considering  all  options  in  the  development  of  a  draft  Onshore 
Order  5. 

Factors  potentially  affecting  gas  variability 

In  addition  to  the  simulation-based  analysis  described  above,  numerous  analyses  were  conducted  on  the 
data  provided  in  an  attempt  to  identify  factors  that  lead  to  a  higher  degree  of  variability.  These  analyses 
included: 

•  heating  value; 

•  temperature; 

•  time  of  production; 

•  pressure; 

•  reservoir  type 

•  type  of  lift;  and, 

•  the  existence  of  a  separator  prior  to  the  meter. 

Heating  value 

For  each  meter,  the  average  heating  value  was  determined  and  then  the  meters  were  separated  into 
one  of  3  categories  depending  on  the  average  heating  value.  The  categories  were  defined  in  order  to 
have  about  the  same  number  of  meters  in  each,  as  follows: 

Low:  less  than  1080  Btu/scf 

Medium:  1081  to  1141  Btu/scf 
High:  more  than  1141  Btu/scf 

For  each  meter,  the  95%-confidence  variation  using  the  student-t  distribution  was  plotted  against  the 
average  time  between  samples.  Meters  with  fewer  than  5  samples  were  excluded  from  the  analysis. 
Figure  9  shows  the  results  of  this  analysis.  Although  it  appears  that  the  medium  category  shows  slightly 
less  variability  overall,  there  is  no  apparent  trend  in  variability  versus  sampling  frequency  for  any  of  the 
three  categories. 

Temperature 

For  this  analysis,  the  sample  temperature  was  plotted  against  heating  value  for  each  meter.  In  order  to 
plot  all  the  data  on  a  single  graph,  both  the  heating  values  and  sample  temperatures  were  normalized  to 
1200  Btu/scf  and  60°F,  respectively.  Because  richer  gases  (higher  heating  value)  are  more  prone  to 
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temperature  effects,  only  the  "High"  heating  value  category  was  used.  Meters  with  less  than  5  samples 
were  excluded  from  this  analysis.  In  addition,  only  39%  of  the  samples  taken  included  sample 
temperature,  thereby  limiting  the  amount  of  data  to  work  with. 

Figure  10  shows  the  plot  of  normalized  heating  value  as  a  function  of  normalized  temperature.  There  is 
no  apparent  dependency  of  temperature  on  heating  value. 

Time  of  production 

The  heating  values  from  meters  in  each  formation  were  plotted  against  the  sample  date  in  order  to 
indentify  common  trends  in  the  data  -  for  example,  if  there  were  rapid  changes  in  heating  value  early  in 
the  life  of  a  well.  Figures  11a  -  lid  show  sample  plots  from  the  same  formations  used  in  Figures  4a  -  4d. 
While  some  meters  exhibit  a  trend  overtime,  there  is  no  commonality  in  the  data.  In  other  words, 
taking  more  frequent  samples  early  in  the  life  of  a  well  does  not  seem  to  be  justified  by  these  graphs. 

Pressure 

For  each  meter,  the  average  static  pressure  was  used  to  categorize  the  meter  as  follows: 

Low:  less  than  50  psig 

Medium:  between  51  and  500  psig 

High:  greater  than  500  psig 

For  each  meter,  the  95%-confidence  variation  using  the  student-t  distribution  was  plotted  against  the 
average  time  between  samples.  Meters  with  fewer  than  5  samples  were  excluded  from  the  analysis. 
Figure  12  shows  the  results  of  this  analysis.  There  is  no  apparent  trend  or  dependency  of  heating  value 
variability  on  static  pressure. 

Reservoir  type 

The  meters  were  categorized  according  to  the  type  of  reservoir  as  follows:  coal  bed  methane,  gas  cap, 
retrograde,  and  tight  sands.  For  each  meter,  the  95%-confidence  variation  using  the  student-t 
distribution  was  plotted  against  the  average  time  between  samples.  Meters  with  fewer  than  5  samples 
were  excluded  from  the  analysis.  Figure  13  shows  the  results  of  this  analysis.  There  is  no  apparent  trend 
or  dependency  of  heating  value  variability  on  reservoir  type. 

Type  of  lift 

The  meters  were  categorized  as  to  the  type  of  lift:  free  flowing,  plunger  lift,  and  pumping  unit.  For  each 
meter,  the  95%-confidence  variation  using  the  student-t  distribution  was  plotted  against  the  average 
time  between  samples.  Meters  with  fewer  than  5  samples  were  excluded  from  the  analysis.  Figure  14 
shows  the  results  of  this  analysis.  There  is  no  apparent  trend  or  dependency  of  heating  value  variability 
on  lift  type. 

Existence  of  a  separator  prior  to  the  meter 

The  meters  were  categorized  as  to  whether  or  not  there  was  a  separator  prior  to  the  meter.  For  each 
meter,  the  95%-confidence  variation  using  the  student-t  distribution  was  plotted  against  the  average 
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time  between  samples.  Meters  with  fewer  than  5  samples  were  excluded  from  the  analysis.  Figure  15 
shows  the  results  of  this  analysis.  There  is  no  apparent  trend  or  dependency  of  heating  value  variability 
on  the  existence  of  a  separator. 

Flow  rate 

The  meters  were  categorized  by  flow  rate  as  follows: 

Low:  Less  than  50  Mcf/day 
Medium:  Between  50  and  100  Mcf/day 
High:  Between  100  and  500  Mcf/day 
Very  High:  Greater  than  500  Mcf/day 

For  each  meter,  the  95%-confidence  variation  using  the  student-t  distribution  was  plotted  against  the 
average  time  between  samples.  Meters  with  fewer  than  5  samples  were  excluded  from  the  analysis. 
Figure  17  shows  the  results  of  this  analysis.  While  the  "Very  High"  category  appears  to  have  less 
variation,  there  were  only  7  meters  in  this  category,  which  is  not  enough  to  make  a  statistically 
significant  conclusion. 
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Variability  vs  Heating  Value 
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Figure  10 


Frontier  Formation 


Fruitland  Formation 

izoo 

nocr 

K  X  x  „  XX  X  xxxxx* 

1000' 

_ T _ 

Figure  lib 

Mesaverde  Formation 

14ULT 

1300" 

*  .  *  a  i  r »  » 

1200" 

x  ff?****  , 

1100" 

■  *  *  ■  ■  ■  ■  m 

1000" 

|  §  1  1  |  1 

|  $  $  $  $  $ 

Figure  lid 


20 


Variability  vs  Pressure 


Variability  vs  Reservoir  Type 


Average  sample  period,  days 
Figure  12 


Average  sample  period,  days 
Figure  13 
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Appendix  1  -  Sample  Calculation  of  Heating  Value  Variability  and  Uncertainty 


To  demonstrate  the  methodology  for  determining  heating  value  variability  (V95%),  a  meter  "S  25-34" 
measuring  gas  from  the  Almy  Formation  will  be  used.  Gas  analysis  data  for  this  meter  was  obtained  for 
the  past  6  years  as  shown  in  columns  (1)  and  (2)  from  Table  1-1  below.  The  average  heating  value  is 
1162.233  Btu/scf  for  the  gas  analyses  shown.  The  standard  deviation  of  the  data  is  calculated  with  the 
following  equation: 


(Eq.  1-1) 

Where 


o  is  the  standard  deviation,  Btu/scf 
HVi  indicates  an  individual  heating  value 
n  is  the  number  of  heating  values  (12  in  this  example) 
is  the  average  heating  value,  Btu/scf 


(1) 

Sample  Date 

(2) 

Heating  Value 
(Btu/scf) 

(3) 

6/26/03 

1153.135 

101.9748 

1/15/04 

1165.056 

3.323116 

6/25/04 

1142.626 

424.6838 

12/11/04 

1154.601 

74.52797 

6/14/05 

1150.726 

156.439 

1/26/06 

1166.713 

12.10373 

5/8/06 

1159.017 

17.77768 

12/4/06 

1179.779 

273.7583 

6/11/07 

1159.117 

16.94358 

12/7/07 

1186.447 

538.8499 

8/29/08 

1179.102 

251.817 

1/15/09 

1162.483 

0.562888 

Average: 

1163.233 

Sum: 

1872.762 

Table  1-1 


For  this  example  the  sum  of  1872.762  is  the  sum  of  the  square  of  the  difference  between  the  individual 
heating  values  and  the  average  heating  value  (Table  1-1,  column  (3)),  so  the  standard  deviation  is: 

=  ±13.05  Btu/scf 
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For  12  samples,  the  variable  tdlst  is  2.201  (see  Appendix  2),  therefore  the  uncertainty  at  a  95%  confidence 
level  is: 


V95%  =  13.05  x  2.201  =  ±28.73  Btu/scf 
Finally,  converting  this  into  a  percentage  gives  a  variability  of: 


=  +2.47% 


From  the  variability,  the  uncertainty  in  average  heating  value  (UHV)  can  then  be  calculated  as  follows: 


Table  1-2  shows  the  uncertainty  in  1-year  average  heating  value  for  30,  60, 180,  and  365  day  sampling 
periods. 


Sampling  Period 

U„v 

(days) 

(±%) 

30 

0.668 

60 

0.944 

90 

1.157 

180 

1.636 

365 

2.330 

Table  1-2 
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2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 


Appendix  2  -  Student  t  distribution  function 


t-dist 

#Samples 

t-dist 

#Samples 

t-dist 

40 

2.023 

160 

1.975 

41 

2.021 

165 

1.975 

12.706 

42 

2.020 

170 

1.974 

4.303 

43 

2.018 

175 

1.974 

3.182 

44 

2.017 

180 

1.973 

2.776 

45 

2.015 

185 

1.973 

2.571 

46 

2.014 

190 

1.973 

2.447 

47 

2.013 

195 

1.972 

2.365 

48 

2.012 

200 

1.972 

2.306 

49 

2.011 

205 

1.972 

2.262 

50 

2.010 

250 

1.970 

2.228 

51 

2.009 

300 

1.968 

2.201 

52 

2.008 

350 

1.967 

2.179 

53 

2.007 

400 

1.966 

2.160 

54 

2.006 

450 

1.965 

2.145 

55 

2.005 

500 

1.965 

2.131 

56 

2.004 

600 

1.964 

2.120 

57 

2.003 

700 

1.963 

2.110 

58 

2.002 

800 

1.963 

2.101 

59 

2.002 

oo 

1.960 

2.093 

60 

2.001 

2.086 

65 

1.998 

2.080 

70 

1.995 

2.074 

75 

1.993 

2.069 

80 

1.990 

2.064 

85 

1.989 

2.060 

90 

1.987 

2.056 

95 

1.986 

2.052 

100 

1.984 

2.048 

105 

1.983 

2.045 

110 

1.982 

2.042 

115 

1.981 

2.040 

120 

1.980 

2.037 

125 

1.979 

2.035 

130 

1.979 

2.032 

135 

1.978 

2.030 

140 

1.977 

2.028 

145 

1.977 

2.026 

150 

1.976 

2.024 

155 

1.975 
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Appendix  3  -  1-Year  Average  Heating  Value  Uncertainty  by  Formation 


Meter  Name 


ANDERSON  COAL 

BAXTER 

BEAR  RIVER/FRONTIER  2 
BIG  GEORGE  COAL 

CANYON  COAL 

CARMEL 

CASTLEGATE 


DAKOTA /MORRISON 
DIETZ  3  COAL 
ERICSON 
FARMINGTON 
FELIX  COAL 
FORT  UNION 
FORT  UNION  COAL 


Meter5 
Meter  6 
Meter  7 
Meter  8 
Meter  9 
Meter  10 
Meter  11 
Meter  12 
Meter  13 
Meter  14 
Meter  15 
Meter  16 
Meter  17 
Meter  18 
Meter  19 
Meter  20 
Meter  21 
Meter  22 
Meter  23 
Meter  24 
Meter  25 
Meter  26 
Meter  27 
Meter  28 
Meter  29 
Meter  30 
Meter  31 
Meter  32 
Meter  33 
Meter  34 
Meter  35 
Meter  36 
Meter  37 
Meter  38 
Meter  39 
Meter  40 
Meter  41 
Meter  42 
Meter43 
Meter44 
Meter45 
Meter  46 
Meter  47 
Meter  48 
Meter  49 
Meter  50 
Meter  51 
Meter  52 
Meter  53 
Meter  54 
Meter  55 
Meter  56 
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Meter  Name 


^samples  Vss% 


FRUITLAND  COAL 


GREEN  RIVER 

HILLIARD/BAXTER 
KNOBLOCH  COAL 
LANCE 


LANCE  POOL 
TRANSILANCE  POOL 
LANCE/WARDELL 


LOWER  CANYON  COAL 
LOWER  COOK  COAL 
MANCOS  B 


MESAVERDE 


MESAVERDE 1 

MESAVERDE  1/LANCE 

MORGAN 

MORRISON 


Meter  57 
Meter  58 
Meter  59 
Meter60 
Meter61 
Meter62 
Meter63 
Meter  64 
Meter65 
Meter  66 
Meter67 
Meter  68 
Meter  69 
Meter  70 
Meter  71 
Meter  72 
Meter  73 
Meter  74 
Meter  75 
Meter  76 
Meter  77 
Meter  78 
Meter  79 
Meter  80 
Meter  81 

Meter  82 
Meter  83 
Meter  84 
Meter  85 
Meter  86 
Meter  87 
Meter  88 
Meter  89 
Meter  90 
Meter91 
Meter  92 
Meter  93 
Meter  94 
Meter95 
Meter96 
Meter  97 
Meter  98 
Meter99 
Meter  100 
Meter  101 
Meter  102 
Meter  103 
Meter  104 
Meter  105 
Meter  106 
Meter  107 
Meter  108 
Meter  109 
Meter  110 
Meter  111 
Meter  112 
Meter  113 
Meter  114 
Meter  115 
Meter  116 
Meter  117 
Meter  118 
Meter  119 
Meter  120 
Meter  121 
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OJO  ALAMO 


PICTURED  CUFFS 


PINON  FRUITLAND  SAND 


POINT  LOOKOUT 


SI-FRONTIER  1/S2-FRONTIER  2 
S1-FRONTIER/S3-BAXTER 
SAN  JOSE 

SHANNON 
TANSILL 
TOCITO 


TRANSIT  ION/MESAVERDE  1 

TURNER 

UINTA 

UPPER  COOK  COAL 


Meter  122 
Meter  123 
Meter  124 
Meter  125 
Meter  126 
Meter  127 
Meter  128 
Meter  129 
Meter  130 
Meter  131 
Meter  132 
Meter  133 
Meter  134 
Meter  135 
Meter  136 
Meter  137 
Meter  138 
Meter  139 
Meter  140 
Meter  141 
Meter  142 
Meter  143 
Meter  144 
Meter  145 
Meter  146 
Meter  147 
Meter  148 
Meter  149 
Meter  150 
Meter  151 
Meter  152 
Meter  153 


Meter  157 
Meter  158 
Meter  159 
Meter  160 
Meter  161 
Meter  162 
Meter  163 
Meter  164 
Meter  165 
Meter  166 
Meter  167 
Meter  168 


WILLIAMS  FORK 
WOLF  CAMP 
WYODAKCOAL 


Meter  169 
Meter  170 
Meter  171 
Meter  172 
Meter  173 
Meter  174 
Meter  175 
Meter  176 
Meter  177 
Meter  178 
Meter  179 
Meter  180 


Sample  Period 


^samples 


V 9 5 %  1-month  3-month  6-month  1-year 


13  5.98 

38  5.53 

13  1.39 

12  1.96 

8  3.18 

12  0.28 

6  5.51 

10  10.25 

19  10.75 

8  0.44 

9  4.91 

7  1.92 

5  1.35 

6  2.13 

11  12.92 

5  23.64 

7  0.59 

8  7.94 

27  13.83 

7  1.71 

6  2.68 

7  0.69 

5  2.01 

7  1.08 

13  1.28 

6  7.85 

5  2.11 

5  1.35 

19  18.73 

18  3.46 

8  7.26 

5  9.68 

11  1.44 

11  0.12 

6  9.72 

21  6.55 

6  6.52 

4  2.15 

5  2.42 

6  1.01 

11  3.38 

11  6.37 

10  0.84 

12  1.40 

12  1.18 

12  1.74 

12  0.81 

12  0.94 

9  1.61 

10  2.77 

4  1.60 

12  6.34 

56  2.81 

42  1.55 

9  0.67 

6  3.44 

6  5.84 

6  0.84 

23  1.08 


1.63  2.83 

1.51  2.61 

0.38  0.66 

0.53  0.93 

0.87  1.50 

0.08  0.13 

1.50  2.60 

2.79  4.84 

2.93  5.07 

0.12  0.21 

1.34  2.32 

0.52  0.91 

0.37  0.64 

0.58  1.00 

3.52  6.10 

6.44  11.16 

0.16  0.28 

2.17  3.75 

3.77  6.53 

0.47  0.81 

0.73  1.26 

0.19  0.33 

0.55  0.95 

0.29  0.51 

0.35  0.61 

2.14  3.71 

0.57  1.00 

0.37  0.64 


1.98  3.43 

2.64  4.57 

0.39  0.68 

0.03  0.06 

2.65  4.59 

1.78  3.09 

1.78  3.08 

0.59  1.02 

0.66  1.14 

0.28  0.48 

0.92  1.60 

1.74  3.01 

0-23  0.40 

0.38  0.66 

0.32  0.56 

0.47  0.82 

0.22  0.38 

0.26  0.45 

0.44  0.76 

0.76  1.31 

0.44  0.76 

1.73  2.99 

0.77  1.33 

0.42  0.73 

0.18  0.32 

0.94  1.62 

1.59  2.76 

0.23  0.40 

0.29  0.51 


4.00  5.69 

3.69  5.26 

0.93  1.32 

1.31  1.86 

2.13  3.03 

0.19  027 

3.68  524 

6.84  9.74 

7.18  1022 

029  0.42 

328  4.67 

1.28  1.83 

0.90  1.28 

1.42  2.02 

8.63  1229 

15.78  22.48 

0.39  0.56 

5.31  7.56 

9.24  13.15 

1.14  1.63 

1.79  2.55 

0.46  0.66 

1.34  1.91 

0.72  1.02 

0.86  122 

5.24  7.46 

1.41  2.00 

0.90  128 

12.51  17.81 

2.31  329 

4.85  6.91 

6.46  920 

0.96  1.37 

0.08  0.12 

6.49  9.24 

4.37  622 

4.35  620 

1.44  2.05 

1.62  2.30 

0.68  0.96 

2.26  321 

4.26  6.06 

0.56  0.80 

0.93  1.33 

0.79  1.13 

1.16  1.65 

0.54  0.77 

0.63  0.90 

1.08  1.53 

1.85  2.63 

1.07  1.52 

4.23  6.03 

1.88  2.68 

1.04  1.48 

0.45  0.64 

2.30  3  27 

3.90  5.55 

0.56  0.80 

0.72  1.03 
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Appendix  4  -  1-Year  Average  Relative  Density  Uncertainty  by  Formation 


ANDERSON  COAL 

BAXTER 

BEAR  RIVER/FRONTIER  2 
BIG  GEORGE  COAL 

CANYON  COAL 

CARMEL 

CASTLEGATE 


DAKOTA /MORRISON 
DIETZ  3  COAL 
ERICSON 
FARMINGTON 
FELIX  COAL 
FORT  UNION 
FORT  UNION  COAL 


Meter  1 
Meter  2 
Meter  3 
Meter4 
Meter5 
Meter6 
Meter  7 
Meter  8 
Meter  9 
Meter  10 
Meter  11 
Meter  12 
Meter  13 
Meter  14 
Meter  15 
Meter  16 
Meter  17 
Meter  18 
Meter  19 
Meter  20 
Meter21 
Meter  22 
Meter  23 
Meter  24 
Meter  25 
Meter  26 
Meter  27 
Meter  28 
Meter  29 
Meter  30 
Meter  31 
Meter  32 
Meter  33 
Meter  34 
Meter  35 
Meter  36 
Meter  37 
Meter  38 
Meter  39 
Meter  40 
Meter41 
Meter  42 
Meter  43 
Meter  44 
Meter  45 
Meter  46 
Meter47 
Meter  48 
Meter  49 
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) 


Formation 


Meter 


^sampler 


Sample  Period 
1 -month  3-month  6-month 


FRUIT  LAND  COAL 

GALLUP 

GREEN  RIVER 

HILLIARD/BAXTER 
KNOBLOCH  COAL 
LANCE 


LANCE  POOL 
TRANSILANCE  POOL 
LANCE/WARDELL 


LOWER  CANYON  COAL 
LOWER  COOK  COAL 
MANCOS  B 


Meter  57 
Meter  58 
Meter  59 
Meter  60 
Meter  61 
Meter  62 
Meter  63 
Meter  64 
Meter  65 
Meter  66 
Meter  67 
Meter  68 
Meter  69 
Meter  70 
Meter  71 
Meter  72 
Meter  73 
Meter  74 
Meter  75 
Meter  76 
Meter  77 
Meter  78 
Meter  79 
Meter  80 
Meter  81 

Meter  82 
Meter  83 
Meter  84 
Meter  85 
Meter  86 
Meter  87 
Meter  88 
Meter  89 
Meter  90 
Meter  91 
Meter  92 
Meter  93 
Meter  94 
Meter  95 
Meter  96 
Meter  97 
Meter  98 
Meter  99 
Meter  100 
Meter  101 
Meter  102 
Meter  103 
Meter  104 
Meter  105 
Meter  106 
Meter  107 
Meter  108 
Meter  109 
Meter  110 

Meter  112 
Meter  113 
Meter  114 
Meter  115 
Meter  116 
Meter  117 
Meter  118 
Meter  119 
Meter  120 
Meter  121 


MESAVERDE  1/LANCE 

MORGAN 

MORRISON 


29 


1.03 


0.27 

0.28 


Formation 


Meter  Name 


^samples  Vs 


Sample  Period 

1-month  3-month  6-month  1-yeai 


OJO  ALAMO 


PICTURED  CLIFFS 


PINON  FRUIT LAND  SAND 


POINT  LOOKOUT 


SI-FRONTIER  1/S2-FRONTIER  2 
S1-FRONTIER/S3-BAXTER 
SAN  JOSE 

SHANNON 

TANSILL 

Tocrro 

TRANSITION 

TRANSIT  ION/MESAVERDE  1 

TURNER 

UINTA 

UPPER  COOK  COAL 


Meter  124 
Meter  125 
Meter  126 
Meter  127 
Meter  128 
Meter  129 
Meter  130 
Meter  131 
Meter  132 
Meter  133 
Meter  134 
Meter  135 
Meter  136 
Meter  137 
Meter  138 
Meter  139 
Meter  140 
Meter  141 
Meter  142 
Meter  143 
Meter  144 
Meter  145 
Meter  146 
Meter  147 
Meter  148 
Meter  149 
Meter  150 
Meter  151 
Meter  152 
Meter  153 
Meter  154 


er155 
sr  156 


WILLIAMS  FORK 
WOLF  CAMP 
WYODAKCOAL 


Meter  157 
Meter  158 
Meter  159 
Meter  160 
Meter  161 
Meter  162 
Meter  163 
Meter  164 
Meter  165 
Meter  166 
Meter  167 
Meter  168 
Meter  169 
Meter  170 
Meter  171 


30 


